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A long-standing puzzle is why there is a difference between the optimal doping ^optimal ~ 0.15 
for the maximal superconducting (SC) transition temperature Tc and the critical doping cScriticai ~ 
0.19 for the highest superfluid density ps in cuprate superconductors? This puzzle is calling for 
an explanation. Within the kinetic energy driven SC mechanism, it is shown that except the 
quasiparticle coherence, pa is dominated by the bare pair gap, while Tc is set by the effective pair 
gap. By calculation of the ratio of the effective and the bare pair gaps, it is shown that the coupling 
strength decreases with increasing doping. This doping dependence of the coupling strength induces 
a shift from the critical doping for the maximal value of the bare pair gap parameter to the optimal 
doping for the maximal value of the effective pair gap parameter, which leads to a difference between 
the optimal doping for the maximal Tc and the critical doping for the highest ps- 
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The parent compounds of cuprate superconductors are 
Mott insulators with an antiferromagnetic long-range or- 
der (AFLRO)^. However, this AFLRO is suppressed by 
doped charge carriers, then superconductivity arises from 
the binding of charge carriers into Cooper pairs'^, thereby 
forming a superfluid with a superconducting (SC) energy 
gap A(k) in the single-particle excitation spectrum. This 
energy gap is corresponding to the energy for breaking 
a Cooper pair of the charge carriers and creating two 
excited states-, while the superfluid density ps is propor- 
tional to the squared amplitude of the macroscopic wave 
fimction^, and therefore describes the SC charge carriers. 
In this case, both A(k) and ps are thus two fundamen- 
tal parameters whose variation as a function of doping 
and temperature provides important information crucial 
to understanding the details of the SC state'^"^. 

After intensive investigations over more than two 
decades, some essential features of the evolution of the 
SC state in cuprate superconductors with doping have 
been experimentally established^^ '^-i where the measured 
energy gap parameter A and the SC transition temper- 
ature Tc show a domelike shape doping dependence, i.e., 
the maximal A and Tc occur around the optimal doping 
i^optimai ~ 0.15, and then decrease in both the under- 
doped and the overdoped regimes^"—. Moreover, the ex- 
perimental measurements^rii throughout the SC dome 
show that the superfluid density ps appears from the 
starting point of the SC dome, and then increases with 
increasing doping in the lower doped regime. However, 
this Ps reaches its highest value around the critical doping 
i^criticai ~ 0.19, and then decreases at the higher doped 
regime, eventually disappearing together with A at the 
end of the SC dome. In particular, it has been showii^"— 
that the maximal Tc around the optimal doping and the 
peak of Ps around the critical doping is a common fea- 
ture of cuprate superconductors. Since A measures the 
strength of the binding of charge carriers into Cooper 
pair*?, while ps is a measure of the phase stiffness^, there- 
fore A and Ps separately describe the different aspects of 
the same SC charge carriers. In this case, a long-standing 



puzzle is why there is a difference between the optimal 
doping for the maximal Tc and the critical doping for the 
highest Ps? 

In this paper, we try to answer this question. Ex- 
perimentally, the measured energy gap A(k) is an ef- 
fective energy gap^^*, which incorporates both the cou- 
pling strength and the bare energy gap A(k). Theoreti- 
cally, the kinetic energy driven SC mechanism has been 
developed^, where Tc is controlled by both the effective 
charge carrier pair gap and the quasiparticle coherence. 
Within this kinetic energy driven SC mechanism, we cal- 
culate the doping dependence of the coupling strength 
Vcs, and the result shows that VcS smoothly decreases 
upon increasing doping from a strong-coupling case in 
the underdoped regime to a weak-coupling side in the 
overdoped regime. Our results also show that the max- 
imal value of the bare charge carrier pair gap parameter 
appears around the critical doping (5criticai ~ 0.195, then 
as a natural consequence, the highest ps occurs around 
this same critical doping. However, the special doping 
dependence of Vos shifts this critical doping for the max- 
imal value of the bare charge carrier pair gap parameter 
to the optimal doping (5optimai ~ 0.15 for the maximal 
value of the effective charge carrier pair gap parameter, 
which leads to that Tc exhibits a maximum around the 
optimal doping. 

Cuprate superconductors have a layered structure con- 
sisting of the two-dimensional Cu02 planes separated by 
insulating layers^. The single common feature is the pres- 
ence of the Cu02 plane, and it seems evident that the 
unusual behaviors of cuprate superconductors are dom- 
inated by this Cu02 plane^.. In this case, it has been 
argued that the essential physics of the doped Cu02 
plane is properly accounted by the t-J model on a square 
lattice^"^. However, for discussions of the difference be- 
tween the optimal doping for the maximal Tc and the 
critical doping for the highest ps, the t-J model can be 
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extended by including the exponential Peierls factors as, 

li)(7 Ifj'a 
la li) 

supplemented by an important on-site local constraint 
'}2a^lcT^i<^ < 1 to remove the double occupancy, where 
the summation is over all sites and for each over its 
nearest-neighbors (NN) 77 or the next nearest-neighbors 
(NNN) 77', C/^ and C;^ are electron operators that re- 
spectively create and annihilate electrons with spin cr, 
S/ = {Sf,Sj,Sf) are spin operators, and ^ is the chem- 
ical potential. The exponential Peierls factors Pj,-, = 
^-i(e/h)A{i)-f, g^jjj p^., ^ ^-,(e/h)A{iyf,' account for the 

coupling of electrons to an external magnetic field in 
terms of the vector potential A(Z)^. To incorporate 
the electron single occupancy local constraint in the t- 
J model ([!]), the charge-spin separation (CSS) fermion- 
spin theoryi^ii^ has been proposed, where a spin-up 
annihilation (spin-down annihilation) operator for the 
physical electron is given by a composite operator as 
C;t ^ ^l-t^r i^U ^ ^li^^)' with the spinful fermion op- 
erator hia = e~**'''/i; that describes the charge degree of 
freedom of the electron together with some effects of spin 
configuration rearrangements due to the presence of the 
doped hole itself (charge carrier), while the spin operator 
Si represents the spin degree of freedom of the electron, 
then the electron single occupancy local constraint is sat- 
isfied in analytical calculations. In this CSS fermion-spin 
representation, the t-J model ll} can be rewritten as, 

ifi' 

- f-i'^hj^hia + Jcs'^Si ■ Si+ij, (2) 

1(7 Ifj 

where Jeff = (1 — ^YJ^ and 5 = (hj^hic) — {h\hi) is the 
doping concentration. 

Since the experimental measurementsii have shown 
that in the real space the gap function and the pairing 
force have a range of one lattice spacing, the hare energy 
gap parameter can be expressed as^^ A = (C'/^C/+f,^ — 

Cl^Cl^f^^) = {hijhi+i^iS+S--hiihi+f,^S-S+^f^). In the 
doped regime without AFLRO, the spin correlation func- 
tions {Si^Si^_^-) ~ {Si^ S'^-) = xi: and then the bare en- 
ergy gap parameter can be rewritten as A = — xiAh, 
with the bare charge carrier pair gap parameter Ah = 
{hi+filhif — hi+fifhii), which shows that the bare energy 
gap is closely related to the bare charge carrier pair gap, 
therefore the essential physics in the SC state is domi- 
nated by the corresponding one in the charge carrier pair- 
ing state. For a microscopic description of the SC state 



in cuprate superconductors, the kinetic energy driven SC 
mechanism has been developecU^ based on the t-J model 
([2]) , where the charge carrier interaction directly from the 
kinetic energy by exchanging spin excitations induces a 
d-wave charge carrier pairing state, and then their con- 
densation reveals the SC ground-state. Moreover, this 
SC state is controlled by both the effective energy gap 
and the quasiparticle coherence. Within this kinetic en- 
ergy driven SC mechanism, the full charge carrier Green's 
function in the zero magnetic field case has been obtained 
explicitly in the Nambu representation asi^^, 

■ \ 7 ^'^nTo +CkT3 - Ahz(k)ri 

where tq is the unit matrix, t\ and T3 are Pauli matri- 
ces, the renormalized charge carrier excitation spectrum 
Ck = ^hF^kj with the mean- field charge carrier excitation 
spectrum = Ztxijk — Zt'x2l\^ — M, the spin correla- 
tion function X2 = 7k = (1/^) Tf, ^"^'^ 7^ = 

{l/Z) Y,fi' , Z is the number of NN or NNN sites, the 
renormalized charge carrier d-wave pair gap Ahz(k) = 
ZhFAh(k), and the charge carrier quasiparticle spec- 
trum i?hk = \l^k + l^hz(k)p, where the effective charge 
carrier d-wave pair gap Ah(k) = Ah(cosfc^ — cosfcj,)/2, 
and is closely related to the self-energy ^,^2^ i\<i, uj) in 

the particle-particle channel as Ah(k) — Ej'^^k, w)|c^=o, 
while the quasiparticle coherent weight ZhF is directly 
associated with the self-energy t!"^^ {k, uj) in the particle- 
hole channel as Zj^p = 1 — ReE£-'(k, = 0)|k=[7r.o]j 
with sio^(k, Lo) is the antisymmetric part of I]^''''(k, w), 

where the self-energies (k, lo) and E2'^-'(k, w) have 
been given in Refs. [H and [H. In this case, the ef- 
fective charge carrier pair gap parameter Ah, ^hF, and 
all the other order parameters have been determined by 
the self-consistent calculatio n^^i^^ . For a convenience in 
the following discussions, the self-consistently calculated 
result^ of Ah versus doping for temperature T = 0.002 J 
with parameters t/J ~ 2.5 and t' /t — 0.3 is replotted in 
Fig. [TJ where the maximal Ah occurs around the opti- 
mal doping ^optimal ~ 0.15, and then decreases in both 
the underdoped and the overdoped regimes. 

With the help of the Green's function ([3]), the bare 
charge carrier pair gap parameter Ah can be evaluated 
explicitly as. 

Ah = ^ ^[cosfc^ - cosfc^]^ ^''^^''^ tanh[^/?j;hk] ■ (4) 

k 

Since the pairing force and this Ah have been incorpo- 
rated into the effective charge carrier pair gap parameter 
Ah^, the strength VcS of the attractive interaction me- 
diated by spin excitations in the kinetic energy driven 
SC mechanism can therefore be obtained in terms of the 
ratio of Ah and Ah as, 

Kff = ^. (5) 
Ah 
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is shown that Tc increases with increasing doping in the 
underdoped regime, and exhibits a maximum aromrd the 
optimal doping, then decreases with increasing doping 
in the overdoped regime, in good agreement with the 
experimental results of cuprate superconductors^^— . In 
particular, Tc that is set by the effective pair gap and the 
quasiparticle coherence has been observed experimentally 
in cuprate superconductors^.. We believe that this prop- 
erty may be a common feature for all superconductors, 
since in spite of the electron-phonon SC mechanism, Tc in 
the conventional superconductors is also determined by 
the effective pair gap and the quasiparticle coherence^S. 



FIG. 1: The effective charge carrier pair gap parameter as a 
function of doping for temperature T — 0.002 J with parame- 
ters t/J = 2.5 and t' /t = 0.3. 
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FIG. 2: (a) The coupling strength and (b) the bare charge 
carrier pair gap parameter as a function of doping for temper- 
ature T = 0.002J with parameters t/J = 2.5 and t'/t = 0.3. 



FIG. 3: The superconducting transition temperature as a 
function of doping with parameters t/J = 2.5, t'/t — 0.3, and 
J = lOOOK. Inset: the corresponding experimental results of 
cuprate superconductors taken from Ref. @. 



Although both Ah and Ah measure the strength of the 
binding of charge carriers into charge carrier pairs. Ah is 
an experimentally measurable quantity, while Ah is not. 
In this case, we have calculated the doping dependence 
of VcB and Ah, and the results of (a) VeS and (b) Ah 
as a function of doping for T = 0.002 J with t/J — 2.5 
and t'/t — 0.3 are plotted in Fig. [51 where the coupling 
strength VoS smoothly decreases upon increasing doping 
from a strong-coupling case in the underdoped regime 
to a weak-coupling side in the overdoped regime, which 
is consistent with the experimental result of cuprate 
superconductorsi^. However, Ah increases with increas- 
ing doping in the lower doped regime, and reaches a 
maximum around the critical doping (^critical ~ 0.195, 
then decreases with increasing doping in the higher doped 
regime. In comparison with the corresponding result of 
Ah in Fig. [H we therefore find that the special doping 
dependence of the coupling strength Vcs in Fig. ^a.) 
induces an important shift from the critical doping for 
the maximal Ah to the optimal doping for the maximal 
Ah, then the doping dependence of Tc is determined by 
Ah (then both Ah and Vcs) and the quasiparticle co- 
herent weight ZhF within the kinetic energy driven SC 
mechanism^^. To see this point clearly, Tc as a function 
of doping with t/J = 2.5, t'/t = 0.3, and J ^ lOOOK 
is plotted in Fig. [3] in comparison with the correspond- 
ing experimental results^ of cuprate superconductors. It 



The essential physics of the domelike shape doping de- 
pendence of Tc in cuprate superconductors can be at- 
tributed to a competition between the kinetic energy (St) 
and magnetic energy (J)^. The parent compounds of 
cuprate superconductors are the Mott insulators^, when 
charge carriers are doped into a Mott insulator, there is a 
gain in the kinetic energy per charge carrier proportional 
to t due to hopping, however, at the same time, the mag- 
netic energy is decreased, costing an energy of approxi- 
mately J per site. As a consequence, the strength of the 
spin excitation spectrum decreases with increasing dop- 
ing, which leads to a decrease of the coupling strength Vcs 
with increasing doping in the framework of the kinetic 
energy driven SC mechanism. Moreover, in the under- 
doped regime, the coupling strength Ves in Fig. O^a) is 
very strong, this implies that the most doped charge car- 
riers can be bound into the charge carrier pairs, then the 
number of the charge carrier pairs and Tc increase with 
increasing doping. However, in the overdoped regime, 
the coupling strength Vcs is relatively weak. In this case, 
not all doped charge carriers can be bound to form the 
charge carrier pairs by this weakly attractive interaction, 
and therefore the number of the charge carrier pairs and 
Tc decrease with increasing doping. In particular, the 
optimal doping is a balance point, where the number of 
the charge carrier pairs and the coupling strength VcS are 
optimally matched. This is why the Tc in cuprate super- 
conductors exhibits a domelike shape doping dependence. 
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Now we turn to discuss the doping dependence of the 
superfluid density. The external magnetic field B — rot A 
applied to the system usually represents a large pertur- 
bation, but the induced field generated by supercurrents 
cancels the external field over most of the volume of the 
sample. As a consequence, the net field acts only very 
near the surface on a scale of the magnetic field penetra- 
tion depth, and then it can be treated as a weak pertur- 
bation on the system as a whole. In this case, the Meiss- 
ner effect can be successfully studied within the linear 
response approach^, where the response current density 
and the vector potential A,^ are related by a nonlocal 
kernel of the response function Kfj^^, as, 

Jt,{<l,^) = - Kf,^{ci,uj)A^{ci,u}). (6) 

i/=l,2,3 

This kernel of the response function in Eq. ([6]) can be 
separated into two parts as Kf^i,{q,uj) = kIu){(i,uj) + 

KjS){q,Lu), where K^^^ and KjS) are the corresponding 
diamagnetic and paramagnetic parts, respectively, and 
are closely related to the current-current correlation func- 
tion. The vector potential A has been coupled to elec- 
trons, which are now represented by Ci^ = hj^Sf and 

Cii ~ ill the CSS fermion-spin representation. In 

this case, the electron polarization operator is expressed 

as P = — e ^ 'R.iCj^Cia = e ^ Ri/i- /i^, and then the cur- 
ler ia 
rent operator j in the presence of the vector potential 
is obtained by evaluating the time-derivative of this po- 
larization operator. According to this current operator j, 
the diamagnetic and paramagnetic parts of the response 
kernel have been obtained in the static limit 

KlrJ i<l,0) = ^5]7p(k + q,k)7:(k-t-q,k)[ii(k,q) 
k 

+ i2(k,q)] =4P)(q,0)5^., (7b) 

where the charge carrier particle-hole parameters 
(f>i = {h\^h,+i^„) and 02 = {h\^hi+fi, „) , A^^ = 
— 4e^(Xi0i< — 2x202^') /^^, while the functions ii(k, q, uj) 
and L2(k, q, w) have been given in Ref.^^. In particu- 
lar, we2^ have shown that in the long wavelength limit, 
i.e., |q| ^ 0, if^a^q -J- 0,0) = at T = 0, reflect- 
ing that the long wavelength electromagnetic response is 
determined by the diamagnetic part of the kernel only. 
However, at T = T^, X^y^(q ^0,0) = -(1/Ai), which 
exactly cancels the diamagnetic part of the response ker- 
nel ([7a|) . and then the Meissner effect is obtained for all 
T <Tc. With the help of the response kernel in Eq. (O, 
the magnetic field penetration depth A(T) by taking into 
account the two-dimensional geometry of cuprate super- 
conductors within the specular reflection model has been 



evaluated as2^, 

oo oo 

A(T) - 4 / h,{x)Ax = - I , ^, (8) 

B } TT J AioKyy(qx,0,0)-Hq2' 



then the superfluid density Ps{T) is obtained as Ps{T) = 
\-\T). 
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FIG. 4: The superfluid density as a function of doping for 
temperature T — 0.002J with parameters t/J = 2.5, t' /t = 
0.3, and J = lOOOK. Inset: the corresponding experimental 
results of cuprate superconductors taken from Ref. [Tol . 

In this case, for the discussions of the difference be- 
tween the optimal doping for the maximal Tc and the crit- 
ical doping for the highest p^, the result of ps as a function 
of doping at T = 0.002 J with t/J = 2.5, t'/t = 0.3, and 
J = lOOOK is replotted in Fig. |3]in comparison with the 
corresponding experimental dataM of cuprate supercon- 
ductors (inset). The result in Fig. |4] shows clearly that 
Ps increases with increasing doping in the lower doped 
regime, and reaches a maximum around the critical dop- 
ing ^critical ~ 0.195, then decreases in the higher doped 
regime, in good agreement with the experimental results 
of cuprate superconductors^""— . In particular, this an- 
ticipated value of the critical doping (^critical ~ 0.195 
is very close to the critical doping ^critical ~ 0.19 ob- 
tained experimentally for different families of cuprate 
superconductors^rJi . 

However, this critical doping (^critical ~ 0.195 for the 
highest Ps is different from that for in Fig. [31 
where the maximal Tc appears around the optimal dop- 
ing ^optimal ~ 0.15. This difference can be understood 
within the present theoretical framework. Since ps is 
related to the current-current correlation function, the 
bare charge carrier pair gap parameter Ah, the coupling 
strength VcS, and all the other order parameters are rele- 
vant, i.e., the variation of ps with doping and temperature 
is coupled to the doping and temperature dependence of 
Ah, Veff, and all the other order parameters^. In par- 
ticular, the doping-derivative of ps at the critical doping 
(^critical ~ 0.195 is obtained as (dps/d(5)|5=5^^.^.^^^| = 0. 
Since ps = X^^, i<ips/dd)\s=5„i^i^^i = is equivalent to 
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{<iX/d5)\s=5^,iti^.^i = 0- III t^is case, {dX/d5)\s=s^,iti^.^i 
can be expressed in terms of Eq. ([5]) as, 



dA 
d^" 



5 = 5critical 



- / dq. 

TT 



dKyy{q^, 0,0) 



[MoKyy(qx,0,0) + q2] 



d^ 



= 0, 



(9) 



5=<5critio 



then it is straightforward to obtain from Eq. (O that 



when (dps/d(5)| 



0, (dAh/d,5) 



0, 



which shows that the doping effects from the coupHng 
strength VcS and all the other order parameters upon 
Ps are canceled each other, then both the maximal Ah 
and the highest ps appear at the same critical doping. 
Moreover, both ps and Ah are the bare quantities and 
separately describe the different aspects of the same SC 
charge carriers. In this case, the domelike shape of the 
doping dependence of ps with the highest value appeared 
around the critical doping is a natural consequence of the 
domelike shape of the doping dependence of Ah with the 
maximal value appeared around the same critical doping 
within the kinetic energy driven SC mechanism. In other 
words, except the quasiparticle coherence, ps is domi- 
nated by the bare charge carrier pair gap parameter Ah, 
while Tc is set by the effective charge carrier pair gap pa- 
rameter Ah, this is why there is a difference between the 
optimal doping for the maximal Tc and the critical doping 
for the highest ps in cuprate superconductors. Finally, we 
have noted that ps dominated by the bare energy gap pa- 



rameter in cuprate superconductors has been observed 
from the photoemission experiments^. Since in the SC 
state, the photoemission peak intensity as a function of 
doping scales with ps, then a measurement of the coher- 
ent component in the quasiparticle excitation has been 
suggested as an indirect measure of the bare energy gap 
parameter in cuprate superconductors^'^. 

In conclusion, within the framework of the kinetic en- 
ergy driven SC mechanism, we have discussed the ori- 
gin of the difference between the optimal doping for the 
maximal and the critical doping for the highest ps in 
cuprate superconductors. By calculation of the ratio of 
the effective and bare charge carrier pair gap parameters, 
we have shown that the coupling strength decreases with 
increasing doping. This special doping dependence of the 
coupling strength induces an important shift from the 
critical doping ^critical ~ 0.195 for the maximal value of 
the bare charge carrier pair gap parameter to the optimal 
doping i5optimai ~ 0.15 for the maximal value of the ef- 
fective charge carrier pair gap parameter, which leads to 
a difference between the optimal doping for the maximal 
Tc and the critical doping for the highest ps in cuprate 
superconductors. 
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